Viscoelasticity of 2D liquids quantified in a dusty plasma experiment 
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The viscoelasticity of two-dimensional liquids is quantified in an experiment using a dusty plasma. 
An experimental method is demonstrated for measuring the wavenumber-dependent viscosity, ri{k), 
which is a quantitative indicator of viscoelasticity. Using an expression generalized here to include 
friction, ri{k) is computed from the transverse current autocorrelation function (TCAF), which is 
found by tracking random particle motion. The TCAF exhibits an oscillation that is a signature 
of elastic contributions to viscoelasticity. Simulations of a Yukawa liquid are consistent with the 
experiment. 
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Two-dimensional (2D) physical systems include elec- 
trons on a liquid helium surface [3, colloids 0, gran- 
ular fluids [H, and dusty plasmas ^. In experiments 
and simulations, elastic properties, such as transverse 
waves [E, 0|, and transport properties, such as viscosity 
V have been studied. 

Viscoelasticity is a property of materials that exhibit 
both viscous and elastic characteristics [Ifl] . One usually 
thinks of viscous properties for liquids and elastic prop- 
erties for solids, but most materials are viscoelastic and 
exhibit both. These include, for example polymers, hu- 
man tissue, and hot metal [lO|. In general, liquids exhibit 
elastic effects especially at short length or time scales ll| , 
but viscous effects at long length or time scales. 

To quantify viscoelasticity, one often uses the 
frequency-dependent viscosity 77(0;) [12], which tends to- 
ward the static viscosity, 77, as a; — 0. The r/(aj) is 
easily measured in three-dimensional (3D) liquids using 
rheometers and viscometers fl3] , but not in most 2D liq- 
uids. 

Besides rj^uj), the wavenumber-dependent viscosity, 
r]{k), has been used by theorists to quantify the viscoelas- 
tic character [l3l - [l7j . They have recently developed ways 
of computing 77 (fc) from the trajectories of random mo- 
tion of molecules However, until now, there have 
been no experimental measurements in any physical sys- 
tems known to us of i]{k) that exploit this new analysis 
method. One difficulty in using this method in an exper- 
iment is that it requires, as its inputs, the positions Xi 
and velocities Vi of N individual molecules or particles 
as they move about randomly. In this Letter, we will 
use an experimental system, dusty plasma, that allows 
observing these inputs directly. 

Here we further develop a method for computing //(fc), 
generalizing it for multiphase systems like dusty plasma. 
As was originally developed for 3D molecular dynamics 
(MD) simulations, the method begins with computing 
the normalized transverse current autocorrelation func- 
fllil (TCAF), which is defined as 



where jy(k,t) = Xlili (0 6xp[i/cxi(i)] is the trans- 
verse current, with the vector k par allel to the x axis. 
Then, ri{k) can be calculated 3, l3] through r]{k)/p = 
l/(fc^$), where $ is the area under the normalized 
TCAF. This equation can be derived, assuming that 
the viscosity is a valid transport coefficient, either from 
the hydrodynamic Navier-Stokes equation or from lin- 
ear response theory ^17]. Here, we generalize this equa- 
tion using dj{r,t)/dt - {rj/ p)V^i{r,t) + Vfi{r,t) = 0, a 
Navier-Stokes equation that includes an additional fric- 
tional drag force J^/j(r, t) due to a second phase 7]. This 
equation is valid in both 2D and 3D systems. Following 
the method of fr?], we find [lit 



7y(fc)/p = ((l/<I>)-7./)/fc2 



(2) 



Here we will characterize viscoelasticity in an experi- 
ment two ways. First, as a signature of elastic effects, we 
wiU detect oscillations in the TCAF [13, Ell for large k. 
Second, we will measure the diminishment of r]{k) as k 
increases. This diminishment occurs along with a rela- 
tive increase of elastic contributions to viscoelasticity, for 
large k. 

Dusty (complex) plasma, is partially ionized gas con- 
taining micron-size particles of solid matter ^4, 19, .20,]. 
Particles have a charge Q and can be electrically confined 
in a single horizontal layer where they self-organize with 
a structure like a crystalline solid [20|. Coulomb repul- 
sion is shielded with a screening length Ad [U, 22|. The 



tion 



CT{k,t) = {j*{k,o)jy{k,t))/{j;{k,o)jy{k,o)), (1) 



elastic properties of the crystalline solid arise from inter- 
particle repulsion and can be characterized by the phonon 
spectrum for longitudinal and transverse waves |23| . 
which have a frequency close to the nominal 2D dust 
plasma frequency ujpd |24l |. The solid can be melted, to 
form a liquid, by laser manipulation [H, [2I) . 

Dusty plasmas are attractive for experimental quan- 
tification of viscoelastic effects at a microscopic scale. As 
in colloids @j and granular fluids [1], they allow video 
microscopy to track the Xi and Vi of individual parti- 
cles. They also provide both elastic and viscous effects. 
The particles are immersed in a medium that is a rar- 
efied gas that does not overdamp particle motion, unlike 
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colloids ^ with their solvents. 

Dusty plasma experiments, until now, have yielded de- 
scriptive presentations of viscoelasticity [27[ and demon- 
strations of the microscopic motion of particles associated 
with viscoelastic response (28| . In experiments, the static 
viscosity has been measured and estimated from dif- 
fusion observations "2^. However, a quantitative charac- 
terization of viscoelasticity, using rjiuj) or ?7(fc), is lacking 
from the literature. 

A challenge in dusty plasma experiments is that they 
do not allow direct contact of the suspension with a con- 
tainer. Thus, the viscoelastic response cannot be mea- 
sured with a rheometer. We overcome this challenge by 
observing the random particle motion and using Eq. ([2]) 
to compute ri{k). We will do this with experimental data, 
and confirm our interpretation using a simulation. 

Using the apparatus of [30| . a plasma was powered 
by 13.56 MHz, 170 V peak to peak voltages. After the 
8.1 /Ltm diameter microspheres were introduced into the 
plasma (which had an Argon pressure of 14 mTorr), they 
experienced a damping rate of — 2.4 s~^ [sij . 

The particles were suspended in a single layer. They 
self-organized in a triangular lattice 20[. Particle mo- 
tion was essentially 2D, with negligible out-of-plane dis- 
placements. The suspension had a diameter w 52 mm 
and contained > 5400 particles. The lattice constant 
b = 0.67 mm corresponds to a Wigner-Seitz radius [24 1 
a = 0.35 mm. 

Particle tracking was done by imaging from the top. 
For each of four runs, 20 s videos were recorded at 
250 frames/s, providing adequate time resolution for the 
TCAF. The (36.2 x 22.6) mm^ field of view (FOV) in- 
cluded w 2100 particles. We recorded the maximum 
5061 frames per run allowed by the 12-bit Phantom 
v5.2 camera, with a lens that provided a resolution 
of 0.03 mm/pixel. For each video frame j, we com- 
puted dll the position of the ith particle, Xij. To com- 
pute jy{k,t), we used Xij = (ijj-i + ^ij' + ^ij'+i)/3 
and vf J — {ijij+i — yij-i)/2St. This finite-difference 
method reduced errors arising from the high frame rate. 
Examples of particle trajectories from the experiment 
are shown in Fig. 1(a). Next, we computed jy{k,t) and 
smoothed its time series over five frames before calculat- 
ing the TCAF, Eq. (P), and finally ry(fc), Eq. ©. 

Before melting the suspension., we used the phonon- 
spectrum method for a lattice [23| to measure Q/e — 
—6000, kq — a/Xo — 0.5, and ojpd = 30 s~^. After 
melting, we determined T from the mean-square velocity 
fiuctuation ^ yielding F = {Q'^/'iT:eoa)/{kBT) = 68. 

We melted the lattice and maintained a steady ki- 
netic temperature T using laser manipulation [2^ |30| . 
Random kicks were applied by radiation pressure from 
a pair of 532-nm laser beams that were rastered across 
the suspension in a Lissajous pattern with frequencies 
= 48.541 Hz and fy = 30 Hz. This pattern filled a 
rectangle larger than the camera's FOV. Along with the 



desired random motion, the Lissajous heating method 
also produces coherent modes [2^, which had about 8% 
of our total kinetic energy for motion in the y direction, 
similar to [s^l ■ We analyzed half of the FOV, where the 
temperature was uniform within extremes of ±20%. 

For comparison to the experiment, we also performed 
a Langevin MD simulation 33-3^ of a 2D Yukawa liq- 
uid to mimic our experiment. Using periodic boundary 
conditions and 4096 particles, the equation of motion 
Eq. (3) of [1^ was integrated, yielding particle trajec- 
tories. Fig. 1(b). The simulation parameters F = 68, 
kq = 0.5, and Vf/upd = 0.08 match the experimental val- 
ues. To improve statistics, the simulation was run much 
longer, ujpdt — 22 300, than the experiment uipdt — 607. 
To validate our Langevin MD simulation, we also per- 
formed a frictionless MD simulation Q and calculated 
r]{k) as in Eq.(l2]) but with Vf = 0; we found that the 
results for ??(fc) for the two types of simulations agree. In 
addition to computing rjik), we also computed the static 
viscosity 77 using the Green-Kubo relation, Eq. (3) of [1]. 
The latter assumes that the shear-stress autocorrelation 
function decays significantly faster than 1/t, which we 
verified. 

Experimental results for the TCAF, Fig. 2(a), reveal 
elastic properties in the viscoelastic regime for this liq- 
uid. The TCAF computed from Eq. ([T]) exhibits an ini- 
tial decay followed by oscillations around zero 0,111, 
for fc& = 3.26 in Fig. 2(a). Such oscillations typically in- 
dicate that the selected wavenumber corresponds to the 
viscoelastic regime. The TCAF is a time series; we also 
calculate its frequency spectrum, shown in the inset of 
Fig. 2(a). (This frequency spectrum can also be used 
in generating a phonon spectrum [ll|). The spectrum 
features a prominent peak at non-zero frequency. This 
peak is a signature of shear elasticity; it would be absent 
in a viscous regime. To our knowledge, the TCAF time 
series has not previously been reported for dusty plasma 
experiments as an indicator of viscoelasticity. 

Simulation results. Fig. 2(b), exhibit features in the 
TCAF and its spectrum |,14] similar to those in the ex- 
periment. This agreement between experiment and sim- 
ulation lends confidence to our use of the TCAF as a 
quantitative indicator of viscoelasticity in an experimen- 
tal system. 

For wavenumbers much smaller than those shown in 
Fig. 2, i.e., for very long wavelengths, we would expect 
viscous behavior characterized by a simple decay of the 
TCAF with no oscillations. This hydrodynamic regime 
has been well studied in simulations and theory Is?! . Ob- 
serving it requires a sufficiently large system. One of the 
attractions of our physical system is that it allows direct 
observation of motion at an atomistic scale. Thus, we use 
it here to observe the viscoelastic regime (at small wave- 
lengths), not the purely viscous hydrodynamic regime. 

As our chief result, our experimentally measured 
wavenumber-dependent viscosity, rj{k), is presented 
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quantitatively in Fig. 3(a). We observe that ri{k) dimin- 
ishes as k increases. PhysicaUy, this trend indicates that 
dissipative or viscous effects diminish at shorter length 
scales. At these shorter length scales, elasticity has a 
greater effect. 

Since previous experiments are not available for quan- 
titative comparison, we compare our experimental results 
to the Langevin simulation. Fig. 3(b). We note that r]{k) 
exhibits the same downward trend and similar quantita- 
tive values in the experiment and the simulations. For 
both the experiment and simulation, we present results 
for r]{k), computed using Eq. ([2|), for the viscoelastic 
regime, i.e., k > 1/b. For each fc, the infinite time limit 
for the integration of 4> was replaced with tj, the time 
of the first upward zero-crossing of TCAF time series 
(Fig. 2). This integration limit retains both the viscous 
effects at short time and the elastic effects within the first 
negative peak. 

Noise in the experimental results arose from the finite 
amount of current data used to compute the TCAF. To 
verify that this accounts for the scatter in the experimen- 
tal r]{k) in Fig. 3(a), we repeated the simulation with a 
shorter time, matching the experiment not only in du- 
ration but also in particle number. This test shows, in 
Fig. 3(b), that scatter arises from the finiteness of the 
jy{k, t) data to the same extent as in the experiment. In 
both the experiment and in the shorter simulation, a few 
TCAF curves were too noisy to analyze, with a lack of a 
well-defined upward zero-crossing; the corresponding few 
data points are omitted from Fig. 3. 

We fit ri{k) in Fig. 3 to the same empirical Fade ap- 
proximant used originally for MD simulations of 3D liq- 
uids of hard spheres [l^ and water [l^. This approx- 
imant, 'r]{k) oc (1 -|- afc^)~^, apparently has never been 
applied for 2D liquids. We found that this form fits both 
our experimental and simulation data in Fig. 3 as well 
as the scatter allows. However, a simple power law does 
not fit the r]{k) data as well. 

In addition to finding that our ri{k) fits the Fade ap- 
proximant, we also find in Fig. 3(b), that it extrapolates 
as fc — J> to the static viscosity rj [l^. In this test, 
we found rj using the Green-Kubo relation Q with our 
Langevin simulation; and this result, shown as a star in 
Fig. 3(b), agrees with previous simulations that used dif- 
ferent methods [1, 0| . 

In conclusion, we performed an experiment to quantify 
viscoelasticity of 2D liquids using the TCAF and r]{k). 
We did this using measurements of random particle mo- 
tion in a dusty plasma, which is a frictional system. We 
generalized a method of calculating r]{k) by including the 
friction in the Navier-Stokes equation; and we presented 
an experimental demonstration of this method. Our ex- 
perimental results for ri{k) show that it diminishes with 
increasing k that can be modeled as oc (l + a/c^)^^, which 
compares well with simulation results. 
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FIG. 1: (Color online) Particle trajectories in a 2D liquid, 
with color representing time. To illustrate the random parti- 
cle motion, (a) shows « 10% of the spatial region we analyzed, 
for a duration 60 uj~]^ which is about ~ 10% of one movie, i.e., 
one run in the experiment, while (b) is a part of a Langevin 
MD simulation, shown over the same time interval. 



FIG. 2: Transverse current autocorrelation function (TCAF) 
in the 2D liquid computed using Eq. ^ for (a) the experi- 
ment at kb = 3.26, and (b) the Langevin MD simulation at 
kb = 3.28. At short times, the TCAF decays due to vis- 
cous effects, while at longer times (after its first positive zero 
crossing, tj) it oscillates due to elastic effects. The frequency 
spectrum for each TCAF, shown in the insets, reveals a peak 
that is a signature of the elastic contribution to viscoelastic- 
ity. These results are different from the pure monotonic decay 
of TCAF and its spectrum that would be observed in a purely 
viscous regime. (Here, b is the lattice constant measured be- 
fore melting.) 



FIG. 3: (Color online) The wavenumber-dependent viscosity 
r;(fc) of the 2D liquid, computed using Eq. ((2)| for (a) the 
experiment and (b) simulations of two sizes. We find that 
ri{k) diminishes with k, which is a signature of viscoelastic 
effects. The size of the smaller simulation mimics the size of 
the experiment; comparing them reveals that the scatter of 
the experimental data (a) arises from the data size. In (b), the 
Green-Kubo (static) viscosity r/ is indicated by a star symbol. 
Here, the kinematic viscosity ri{k)/p and wavenumber k are 
normalized to be dimensionless. 
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